A Family of One-Dimensional Vlasov- Maxwell Equilibria for the 



A family of self-consistent collisionless distribution functions for the force-free Harris sheet is pre- 
sented. This family includes the distribution function recently found by Harrison and Neukirch 
[Phys. Rev. Lett. 102, 135003 (2009)] as well as distribution functions with a different dependence 
on the particle energy, but with the same dependence on the canonical momenta. It is shown gen- 
erally that the other distribution functions in the family give rise to the same pressure function and 
thus to the same current density and magnetic field as the known distribution function, provided 
certain conditions on the parameters are satisfied. A number of examples of distribution functions 
from the new family are given, which illustrate the use of the general method. 
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I. INTRODUCTION 



Investigations of plasma instabilities and plasma waves, both in astrophysics and in the 
laboratory, frequently start with a consideration of equilibrium solutions. For collisionless 
plasmas, the required equilibria are found by solving the steady-state Vlasov-Maxwell equa- 
tions. Collisionless plasma equilbria have been investigated for many years, and a number 



of examples can be found in Refs. llNlQ 



Force-free magnetic fields, which satisfy j x B = 0, such that the current density and 
magnetic field are parallel to each other, can be used to model magnetic fields in low-beta 
plasmas, such as that of the solar corona. Finding collisionless force-free equilibria is, how- 
ever, not a trivial task, and hence there are few known examples. Of these known examples, 
only one is of the nonlinear force- free type^^J^ (the rest are linear force-free^'^'^'^) . 

It is well known that the steady-state Vlasov equation has many solutions. In fact, any 
positive function depending only on the constants of motion is a solution, provided the ve- 
locity moments exist. It can be shown that Ampere's law, V x B = fiQ^, can be written in 
terms of derivatives of the pressure function with respect to the components of the vector 
potential. An obvious consequence of this result is that two distribution functions giving 
rise to the same pressure function will automatically satisfy the Vlasov-Maxwell equations 
for the same magnetic field profile. 

In the present paper, it is shown how this property can be used to calculate a new family 
of distribution functions for the force-free Harris sheet. This family includes the previously 
known distribution function^^ii^, as well as more general distribution functions with a dif- 
ferent dependence on the particle energy, but with the same dependence on the canonical 
momenta. A direct comparison can be made between the pressure function calculated from 
the more general distribution functions and the pressure function calculated from the known 
distribution function in order to give conditions on the parameters which, when satisfied, 
will mean that the general distribution functions are also valid solutions for the force-free 
Harris sheet. 

The idea of finding new distribution functions for a given magnetic field profile by chang- 
ing the dependence on the particle energy has been used before, but only for cases depending 
on the particle energy and a single component of the canonical momentum. Fu and HaiJ^, 
for example, showed that kappa type distribution functions can be used for the Harris sheet 
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magnetic field profile^. More recently, Kocharovsky et al.— discussed distribution functions 
with an arbitrary dependence on the particle energy and a fixed dependence on one compo- 
nent of the canonical momentum for the relativistic case. It must, however, be emphasized 
that, for finding force-free Vlasov-Maxwell equilibria, it is crucial that the distribution func- 
tions depend on two components of the canonical momentum, and that the magnetic field 
has more than one non-zero component .-'^'^'i^'i^"—. 

In the present paper, the general theory behind the work is given in Section [Til then the 
force-free Harris sheet model is discussed in Section IIIII Section [IV] contains a discussion 
of the method used to find the new equilibria, and this method is illustrated in Section |Vl 
where three examples are given. Finally, Section IVII contains a summary and conclusions. 



II. GENERAL THEORY 



Finding Vlasov-Maxwell equilibria involves solving the steady-state Vlasov equation, 

+ — (E + vxB)-^ = 0, (1) 
together with the steady-state Maxwell equations, 

V ■ E = -, (2) 

V X E = 0, (3) 

V X B = /ioj, (4) 
V-B = 0. (5) 

In the present model, it is assumed that each quantity varies only with the 2;-coordinate, 
so that there is spatial invariance with respect to the x- and y-coordinates, and also time 
independence (since the system is in equilibrium). Furthermore, it is assumed that the 
magnetic field vanishes in the ^;-direction, and that it can be written as B = V x A, 
where A = {A^, Ay, 0) is a vector potential whose ^-component also vanishes. The x- and 
y-components of B are then given by 
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and so Eq. ([5]) is clearly satisfied. It is also assumed that the electric field E can be written 
as E = — V0, where (f) is an electric potential, so that 




and so the electric field clearly satisfies Eq. ([3]). A further assumption is that the plasma 
consists of two particle species (ions and electrons) with charges Qi = e and qe = —e. 

Due to the symmetries of the system, there are three constants of motion, namely the 
Hamiltonian (particle energy), Hg, arising from the time independence of the system, given 
by 

Hs = ^rusivl + vl + vl) + qs(j), 

the canonical momentum in the x-direction, p^g, arising from the spatial symmetry in the 
x-direction, given by 

Pxs = rrigV^ + qsA^, 

and the canonical momentum in the y-direction, pyg, arising from the spatial symmetry in 
the y-direction, given by 

PyS = TTlgVy + qgAy, 

where rUg is the mass of particle species s. 

Solving Eq. ([T]) by the method of characteristics gives a set of functions fg = 
fs{Hg,pxg,Pyg), which depend only on the constants of motion. These are distribution func- 
tions for species s. The remaining Maxwell equations to be solved are Gauss' law (Eq. ([2])) 
and Ampere's law (Eq. (j4])), which can be expressed as 

1 d^A^ 

1 d^Ay _ , 

/io dx^ 

where the charge density a and the x- and y-components of the current density, jx and jy, 
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can be expressed as 

/OO Qp 
fs{Hs,Pxs,Pys)d^V = (8) 

/OO gp 
VxfsiHs,Pxs,Pys)d^V = — (9) 
^ -OO Cl/ix 

/OO Qp 
Vyfs{Hs,Pxs,Pys)d^V = — (10) 

where Pzz is the 2;2;-component of the pressure tensor, defined as 

Pzz{Ax,Ay,<j)) = Y,mJ vlfsd\. (11) 

Throughout this paper, strict charge neutrahty (0 = 0) will be assumed, in line with 
previous work-»i^"— , leading to the condition a = 0. Ampere's Law can then be written 

/"o^, (12) 



'-X 

dz^ dAx 

d^Ay dPzz „S 

-d^ = m;- 

If the distribution functions are known, it is relatively straightforward to calculate the mag- 
netic field profile by solving Eqs. ( fT2|) and ( |T3l) . However, determining a distribution func- 
tion for a given magnetic field profile is much more difficult and involves solving Eq. ( ITTll 
as an integral equation^'^. Distribution functions giving rise to the same pressure func- 
tion Pzz{Ax, Ay) will give rise to the same current density through Eqs. IQ and f|TO|) . and 
hence the same magnetic field profile through Eqs. (fT2!) and (fT3|) . This idea is central to 
the method employed in Section IIVI to find new equilibrium distribution functions for the 
force-free Harris sheet. 



III. THE FORCE-FREE HARRIS SHEET 



The force-free Harris sheet has a magnetic field profile as follows, 

Bxjfhs = Sotanh(z/L), 

R - ^° 

^■^^"^ ~ cosh{z/L)' 

where Bq is a constant and L is a parameter which specifies the thickness of the sheet (as 
mentioned in Section [III)- The x-component of the field is the same as that of the Harris 
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sheetj^, and the addition of the y-component makes the field force-free, since Bl + = Bq. 
The ^/-component of the field also maintains the force balance across the sheet, since both 
the plasma pressure P^z and the magnetic pressure -B^/2/io = {B^ + By) /2^q are constant 
with respect to z^'—, but the magnetic field components vary with z. Fig. [T] shows a plot 
of the field lines for the force- free Harris sheet. Using Eqs. and ([7]) gives the components 




FIG. 1: Field line plot for the force- free Harris sheet, 
of the vector potential as 

A^jfhs = 25oLtan-i(e"/^), 
Ajfhs = - BqL In [cosh{z / L)] . 

The current density can be calculated from Ampere's Law, and is given by 

Bq tanh(2;/L) 

•^"•^^^^ " ;;i^cosh(^/L)' 

_ ^0 1 

^^•^^'^ fioL cosh' {z/LY 
Fig. [2] shows the magnetic field and current density profiles for the force- free Harris sheet. 
The magnetic field is a nonlinear force-free field, and so the current density is parallel to the 
magnetic field, 
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FIG. 2: Plot of the magnetic field and current density profiles for the force- free Harris sheet. 



with the force- free parameter a{z) given by 



a{z) 



L cosh(z/L) 

which is constant along a given field line, but varies from field line to field line. 

A distribution function satisfying the steady-state Vlasov-Maxwell equations for the force- 
free Harris sheet, which was found by Harrison and Neukirch^, is given by, 

nos 



fsjfhs 



exp{-l3sHs) [as cos{l3sUxsPxs) + exp{l3sUysPys) + bs 



(14) 



where as,bs,Uxs,Uys,nos and (3s are constant parameters of the distribution function, and 
Vth,s is the thermal velocity of particle species s. Note that the condition bs > as must hold 
to ensure that fs is positive. The zz-component of the pressure tensor, Pzzjfhs, is given in 
terms of the microscopic notation by 



zzjfh. 



■ cos{ePeUxeAx) + exp{-ePeUyeAy) + b 



19 



(15) 
. In the 



A detailed discussion of the derivation of Eqs. ffT^ and (flSj) can be found in Ref. 
next section, it is shown how to calculate a family of distribution functions for the force-free 
Harris sheet, which includes the distribution function iHM . as well as distribution functions 
with a different dependence on the particle energy Hs- 
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IV. FINDING A FAMILY OF DISTRIBUTION FUNCTIONS FOR THE FORCE- 
FREE HARRIS SHEET 



Consider a more general distribution function of the form 

TTV^ ( Tfl a \ 

= ^/o^ ( \ ^ 1 [«s COs(/3sUxsPxs) + Q-^v{fisUysVys) + &s], (16) 

where h is an arbitrary function of {nisH g) / qj. and /o is a positive constant. The distribution 
function f ll6p is of the same form as the distribution function given by Eq. f ll4p . but it has a 
different dependence on the particle energy. It can be shown that, when calculating velocity 
moments of distribution functions of the form 

it is possible to make the resulting integrals independent of the particle species s whenever 
= 0, that is, whenever strict charge neutrality is assumed. This has been shown before 
for distribution functions depending only on energy and one component of the canonical 
momentum^i^ and can be used to show that distribution functions of the form f|T6l) are 
solutions of the Vlasov-Maxwell equations for the force- free Harris sheet, provided certain 
conditions on the parameters are satisfied. 

The 2;2;-component of the pressure tensor, Pzz,g, can be written as 



where d?p = dp^gdPysdPzs- At this point, the dependence of Pzz,g on is still stated explicitly. 
In the later parts of the paper, however, the assumption = will be made. 
Defining three new variables E, P and Q as 

E = 

p Pxs 

Q = ^, 

and using these as integration variables, gives the pressure as 

I I 1*00 rOO rOQ 

„ J -co J -co J Emin 



2m„0 , , ^^ 

2E ^ - (P - A^f - (g - Ayf 

Is 



1/2 

dEdPdQ, (18) 
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where 



2^. ■ .... 

The charge density, cTg, can be calculated from Pzz,g using the definition in Eq. ([8]), and is 
given by 

1^ I POO roo roo 
Qs J -IX) J -co J E^A„. 



X 



2m 



-OO J —CO J Emir, 

n-1/2 



2 . )2 



2S - (P - A,)' - (Q - 



dEdPdQ. (19) 



It can be seen that the terms inside the integration in Eq. f[T9|) will all be independent of s, 
that is they will not depend upon the particle species, if the following three conditions are 
satisfied, 

= 0, (20) 

e/3e|Mxe| = el3i\u^i\ = a, (21) 
-e/3eUye = e(3iUyi = 7. (22) 

The above conditions are consistent with those obtained in Ref. Il9|. Due to the modulus 
signs in Eq. ( 12T]) . the parameter a is always positive, and so Uxe and Uxi from the distribution 
function f fT6|) can be positive or negative, and can have the same or opposite sign from each 
other. The neutrality condition a = gives 



where 



00 rco roo 



h = / / h{E)cos{aP) 

J —CO J —00 J Emin, 

X [2E - (P - Axf - (Q - Ayf] dEdPdQ, 

/OO POO poo 
/ / h{E)exp{^Q) 
-co J —CO J Emin,0 

X [2E - (P - Ax)^ - (Q - AyY] dEdPdQ, 

/OO POO poo 
/ / HE) 
'OO J —CO J Emin, 

X [2E - (P - Axf -{Q- AyY] '^''^ dEdPdQ, 
where -Emm,o is the value of Emin when (J) = 0. This then gives 



[tti - ae)h + {hi - he)h = 0, 



which will of course be satisified if 



ae = ai = A, (23) 
he = h = B. (24) 

Note that the condition hs > as mentioned in Section IIIII leads to the condition B > A, 
which must be satisfied to ensure that the distribution functions are positive. When the 
conditions (120|) - (!22|) . (123!) and are satisifed, the pressure Pzz.g, given by Eq. (fTSl) . can 
be rewritten as 

(too _|_ irn ■ ) p 

PzzM-^ ^y) = Vo— —[AH, cos(«A,) + H2 eMl^y) + BH^], (25) 

where 

/ / h{E) [2E -S^- T^] ^''^ cos{aS)dEdSdT, 

-OO J — OO J Emin,Q 

/OO poo roo 
/ / h{E) [2E -S^- T^] ^''^ expi-fT)dEdSdT, 
-OO J —OO J Emin,0 

/OO POO POO 
/ / h{E) [2E -S^- T^] ^''^ dEdSdT, 
-OO J —OO J Emin,0 

with S = P — Ax and T = Q — Ay. This can be compared with the pressure (ITSll . and so 
the general pressure given by Eq. ( l25l) will be equal to the pressure (fTSll if the following 
additional conditions are satisfied, 

2(me + mi)e /3e + A 

= no, (26) 

^ (27) 

„ (28) 

When these conditions are satisfied, in addition to conditions f l20|) - fl22|) . fl23|) and flM|) . 
Ampere's law in the form given by Eqs. (fT2|) and (1131) is satisfied for the set of general 
distribution functions f lT6|) and, therefore, they form a family of equilibrium solutions of the 
Vlasov-Maxwell equations for the force-free Harris sheet, in addition to the known distribu- 
tion function^*. In the next section, three explicit examples are given, which show possible 
choices of the function hijUsHs/ q^) = h{E). The validity of the conditions (I26l)-(l28ll will of 
course depend on the choice of the function h{E), and this will be discussed in each of the 
three examples. 
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V. EXAMPLES OF NEW DISTRIBUTION FUNCTIONS 



The following three examples in this section illustrate the use of the method discussed in 
Section [TV] for various choices of the function hijUsHg/ ql) = h{E). 



A. Delta Function 



Consider a distribution function of the form fll6p . with the function h{msHs/ q^) = h{E) 
given by a delta function, 

h{E) = 6{E-Eo), (29) 

where £'0 > Emin,o- This corresponds to a case where the distribution function is zero 
everywhere except for one particular value of the energy {E = Eq), and so all particles are 
assumed to have the same energy. Carrying out the i?-integration first gives Pzz,g as 



/OO noo 
/ (2^0 -S^- T^)^/^ cosiaS)dSdT 
-00 J —CO 

/OO poo 
/ {2Eo -S^- T'^fl'^ exp{^T)dSdT 
-00 J —00 

/OO poo 
/ {2Eo -S^- T'^f'^dSdT 
-00 J —00 

Using cylindrical coordinates (r,^), with 

^2 + 2-2 _ 2Eor^ 



S = \'2EQr cos9, 
T = ^/2E^rsme, 



(30) 
(31) 
(32) 



then gives 



(me + mi)e 

"271 /"I 



/'2n pi 

Acos{aArc) / / r(l — r^)"*^/^ cos(q;V cos 6')(ir(i6' 
Jo Jo 



+ exp{'-^Ay) / / r (1 — r^)"*^/^ exp(7V sin 6*) (ir(i6' + 
Jo Jo 



2'kB 



where 



a 



i = V2^7- 



(33) 
(34) 
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The remaining integrations can then be carried out by using the formulae (lAip -( IX4|) given 
in the Appendix. This gives Pzz,g as 

(me + mi)e 



rrtprnj 



A 



TT 



TXl 



B 



Jsnia') cos(aA^) 



(35) 



The Bessel functions J3/2 in Eq. f l5^ . of fractional order, can be written in terms of spherical 
Bessel functions ji (of integer order) through the identity 



TT 



Jn[Z) 




2z 



(36) 



which gives, 



where 



This gives Pzz,g as 



^3/2(2;) 



2z 




TT 



sm z cos z 



(37) 



Pzz,g{A^,Ay) = 47r/o 



meTrii 7 



/3 



(7' cosh 7' — sinh 7') 



7'^ / sin a' — a' cos a' 
A coiiiaAr)—^ 



a'^ \7' cosh 7' — sinh 7' 
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+ exp(7Ay) + — 



7 



/3 



(38) 



3 \ 7' cosh 7' — sinh 7' ^ 

The conditions (!26l) -( l28l) for the pressure (!38l) to be equal to the pressure (|T5l) are then given 
by 



no 


= 47r/o^ 


1 


- A^" ( 


2 
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3 V7 



"'-^i^v:^ (7' cosh 7' -sinh 7') 



sin a' — a' cos a' 



7' cosh 7' — sinh 7' 

/3 



7' 



(39) 
(40) 
(41) 



and the condition B > A gives rise to the following condition on b 

>/3 



6>i 



6 \ sin a' — a' cos a' 



(42) 
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The right-hand side of condition ( 139|) is always positive since the 7'-dependent function, 
(7' cosh 7' — sinh7')/7'^, is always positive, regardless of the value of 7'. Note also that 
rris, e, Eq, /q > 0. Condition ( 15^ is, therefore, a valid condition since its left-hand side is 
also always positive (/3s, rig > 0). The constant A in condition fHU]) is positive, and the 
7'-dependent part is the same as in condition f l39p which, as discussed, is always positive for 
any value of 7'. The a'-dependent part is given by 



sin a' — a' cos a' ji {a') 



a 
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a' 



(43) 



which can be positive or negative due to the spherical Bessel function ji{a'), shown in Fig. 
m Condition fHOj) is, therefore, only valid in the regions where ji{a') > 0. It can be seen 
from Fig. [3] that it is valid for small values of a'. Condition fHTj) is also a valid condition, 
since both the left- and right-hand sides are positive {B is positive and the 7'-dependent 
part is positive). The fact that 7' can be positive or negative means that the parameters Uye 
and Uyi from the distribution function ( |T6l) (with the function h{msHs/ q^) = h{E) given by 
Eq. fl29l) ) can be positive or negative, but must have opposite signs from each other through 
the definition ([22]) of 7 = ij ^[Wq. 



0,50 



0,25 



0,00 



-0.25 




FIG. 3: The spherical Bessel function ji(a') 
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B. Step Function 



Consider a distribution function of the form (fT6|) . with the function h{msHs/ q^) = h{E) 
given by a step function, 



h{E) = Q{Eo-E) 




(44) 



The ii^-integral is the same for each part of Pzz,g, and is given by 

POO pEo 

/ Q{Eo- E){2E - -T^ f^dE = {2E - - T^f'^dE 

^ E 771171,0 ^ ^min,0 

= -{2Eo-S^-Ty/\ 
3 

Using a cyhndrical coordinate system {r,9) as in the first example (see Eqs. (150]) the 
resulting integrals can then be evaluated in a similar way to the previous example (see the 
Appendix for details). Then, using Eq. fl36|) as in the previous example, the resulting Bessel 
functions (of order 5/2) can be expressed in terms of the spherical Bessel function j2, where 



The pressure is then given by 



— ] smz cos z. 

z-^ z I z^ 



Pzz,g{A^,Ay) = AttJo — ^^^^ — [(3 + 7 )smh7 -87 cosh7j 



meTTli 7" 
,/5 



A cosiaA. 



7 (3 — a ) sin a' — 3a' cos a' 
^ a'^ (3 + 7'^) sinh 7' — 37' cosh 7' 



+ exp(7ylj^) + 



B 



7 



/5 



15 (3 + 7'^) sinh 7' — 37' cosh 7' 



(45) 



where a' and 7' are given by Eqs. ( 133|) and The conditions ( 126|) -( |28|) for the pressure 
[5]) to be equal to f|T5|) are then given by 



„/5 



——no = Anfo 

PePi merrii 7 

1 ^'j'^ (3 — a'^) sin a' — 3a' cos a' 

2 a'^ (3 + 7'^) sinh 7' — 37' cosh 7' 



5 



7' 



/5 



15 (3 + 7'^) sinh 7' — 37' cosh 7' ' 



(46) 
(47) 
(48) 
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and the condition B > A gives the following condition on b, 

b ^ "^'^ 49 

30 (3 — a'^) sin a' — 3a' cos a' 

The right-hand side of condition f l46l) is always positive since the 7'-dependent function, 

((3 + 7'^) sinh7' — 37' cosh 7') /7'^, is always positive, regardless of the value of 7'. Note also 

that rris, e, Eq, /q > 0. Condition is, therefore, a valid condition since its left-hand side 

is also always positive (/3s, > 0). The constant A in condition fH7|) is positive, and the 

7'-dependent part is the same as in condition ( H6|) which, as discussed, is always positive for 

any value of 7'. The a'-dependent part is given by 

(3 — a'^) sin a' — 3a' cos a' J2(q^') 



a'^ a'^ 



(50) 



which can be positive or negative due to the spherical Bessel function j2{c('), which has 
a similar profile (qualitatively) to that of the spherical Bessel function ji{a'), shown in 
Fig. |3l Condition fH7|) is, therefore, only valid in the regions where j2(tt') > 0. Condition 
(HHj) is also a valid condition, since both the left- and right-hand sides are positive {B is 
positive and the 7'-dependent part is positive). The fact that 7' can be positive or negative 
means that the parameters Uye and Uyi from the distribution function 0161) (with the function 
h{msHs/ql) = h{E) given by (HI])) can be positive or negative, but must have opposite signs 



from each other through the definition ( |22l) of 7 = 7'/ \^2Eq. 



C. Power of Eo- E 

Consider a distribution function of the form (fT6|) . with the function h{msHs/ q^) = h{E) 
given by 



h{E) 



{Eo - E)^, E<Eo 
0, Eo<e' 

where x > — 1- When calculating Pzz,g using Eq. ([2SD, the ii^-integration is the same in each 
triple integral, and is given by 

/ {Eo-E)^[2{E-Ernra,o)Y^^dE, (51) 

^ Emin,0 

which can be written as 

V2 I ^^{ilJo-^Y'^diP, (52) 
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where ip = Eq — E and ipo = Eq — Emin,o- The integral ( 152|) can be evaluated by using the 
formula 

[ t''-\z-tY-^dt = z''+''-^B{ix,u), 
Jo 

where -B(/i, z/) is a beta function defined by 

B{^,iy) = — ■ — r, 
r(/i + z/) 

and 3f?/i > 0, 3?z/ > (which gives x > The integral (IHTj) is then given by 

4.2xr + 5/2^ " ' 

(note that Emmfl = (1/2)('S'^ + T^)). Using a cylindrical coordinate system as before (see 
Eqs. (EOD-ia)) gives P,,,g as 

'^^ ' memi 2x+ir(x + 5/2)^ 

r-27r A-l 



Acos{aA.^) / r(l -r^)^/^+^ cos(aV cos 6')rfrrf6' 
-/o Jo 

+ exp{'yAy) / / r(l - r^)^/^+'^ exp(7V sin 6')(ir(i6' H 

Jo Jo X + 5/2 



where, as in the previous two examples, a' and 7' are given by Eqs. f p3|) and fl34|) . As in the 
previous two examples, the remaining integrations can be carried out by using the formulae 
in the Appendix. The pressure Pzz,g is then given by 



^(%y ^^^±^cos(aA.)+exp(7A 

V" / ^5/2(^7') 

5(^y)x+5/2 



2x+3/2r(^ + 7/2)J^^5/2(^7') 



(53) 
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The conditions (126|) -( 128|) for the pressure (153!) to be equal to (fT5|) are then given by, 

PePi meTTli 

2 V«'/ ■/x+5/2(n') ^ ^ 



2x+3/2^(;^ + 7/2)J^+5/2(^7' 
Conditions (15^ to (l56l) reduce to the two conditions 

'>''x+5/2(y 



(56) 



^'X+5/2 



> 0, (57) 



%^ > 0, (58) 

on removing known positive quantities from the inequahties and using the fact that Jy{iz) = 
Jy{z). The condition B > A gives the following condition on h 

(aO^+5/2 

> 2X+5/2r(;^ + 7/2)J,+5/2(«0' 

The condition ( 1571) contains a root of 7' meaning that, in general, 7' must be positive, and so 
the conditions flFTI) and fISS]) are satisfied in the regions where 7^+5/2(0;') and ^^+5/2(7') are 
positive (note again that a' = > through the definition (I2T1) ). Note that, when x 

is an integer, the identity (!36|) can be used to express the Bessel functions Jn+1/2 in terms of 
spherical Bessel functions j„. This was done in the previous example (which corresponds to 
X = 0). In that example, the conditions ( 154|) to ( |56|) were expressed in terms of sin a', cos a', 
sinh7' and cosh 7' (from the spherical Bessel function ji given by ( 1371) ) and it was observed 
that 7' could be positive or negative without violating the conditions on the parameters (the 
fractional part of the power of 7' cancelled out). 



VI. SUMMARY AND CONCLUSIONS 



A method has been presented, which shows how to calculate a new family of distribution 
functions for the force- free Harris sheet, providing certain conditions on the parameters 
are satisfied. This family includes the known distribution function found by Harrison and 
Neukirch^^, as well as distribution functions with a different dependence on the particle 
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energy, but with the same dependence on the canonical momenta. Three specific examples 
have been given of possible solutions, to illustrate how the method can be used. In each 
example, the conditions on the parameters have been stated explicitly, and the validity of 
these conditions has also been discussed. So far, the force-free Harris sheet is the only 
magnetic field profile for which nonlinear force-free Vlasov-Maxwell equilibria are known. 
If an equilibrium solution was discovered for another field profile, then this method could 
potentially be used to extend the known solution to a whole family of solutions for that 
particular field profile. 
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APPENDIX: EVALUATION OF INTEGRALS 

In each of the three examples in Section |Vl after changing to a cylindrical coordinate system 
(r, 9), the ^-integrations can be carried out by using the formulae 



where Jq and /q are Bessel functions of the first and second kind, respectively. The r-integrals 
then take the following form. 



^0 

Integrals of the form ( ]A3[) can be evaluated firstly by using the substitution r = sint, and 
then by using the formula 




(Al) 



(A2) 





which is valid for 3fJ/x, 3fJz/ > — 1. This gives 




(A3) 
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Integrals of the form (1A3|) can be evaluated by firstly using the identity 

/o(a) = Jo{ia), 

and then by using the same steps that were used to evaluate the integrals of the form flASp . 
This gives 
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